A note on the use of the generalized odds ratio in meta-analysis of association studies involving bi- and tri-allelic polymorphisms

Abstract Background The generalized odds ratio (GOR) was recently suggested as a genetic model-free measure for association studies. However, its properties were not extensively investigated. We used Monte Carlo simulations to investigate type-I error rates, power and bias in both effect size and between-study variance estimates of meta-analyses using the GOR as a summary effect, and compared these results to those obtained by usual approaches of model specification. We further applied the GOR in a real meta-analysis of three genome-wide association studies in Alzheimer's disease. Findings For bi-allelic polymorphisms, the GOR performs virtually identical to a standard multiplicative model of analysis (e.g. per-allele odds ratio) for variants acting multiplicatively, but augments slightly the power to detect variants with a dominant mode of action, while reducing the probability to detect recessive variants. Although there were differences among the GOR and usual approaches in terms of bias and type-I error rates, both simulation- and real data-based results provided little indication that these differences will be substantial in practice for meta-analyses involving bi-allelic polymorphisms. However, the use of the GOR may be slightly more powerful for the synthesis of data from tri-allelic variants, particularly when susceptibility alleles are less common in the populations (≤10%). This gain in power may depend on knowledge of the direction of the effects. Conclusions For the synthesis of data from bi-allelic variants, the GOR may be regarded as a multiplicative-like model of analysis. The use of the GOR may be slightly more powerful in the tri-allelic case, particularly when susceptibility alleles are less common in the populations.

Functional markers for gene mapping and genetic diversity studies in sugarcane

Abstract Background The database of sugarcane expressed sequence tags (EST) offers a great opportunity for developing molecular markers that are directly associated with important agronomic traits. The development of new EST-SSR markers represents an important tool for genetic analysis. In sugarcane breeding programs, functional markers can be used to accelerate the process and select important agronomic traits, especially in the mapping of quantitative traits loci (QTL) and plant resistant pathogens or qualitative resistance loci (QRL). The aim of this work was to develop new simple sequence repeat (SSR) markers in sugarcane using the sugarcane expressed sequence tag (SUCEST database). Findings A total of 365 EST-SSR molecular markers with trinucleotide motifs were developed and evaluated in a collection of 18 genotypes of sugarcane (15 varieties and 3 species). In total, 287 of the EST-SSRs markers amplified fragments of the expected size and were polymorphic in the analyzed sugarcane varieties. The number of alleles ranged from 2-18, with an average of 6 alleles per locus, while polymorphism information content values ranged from 0.21-0.92, with an average of 0.69. The discrimination power was high for the majority of the EST-SSRs, with an average value of 0.80. Among the markers characterized in this study some have particular interest, those that are related to bacterial defense responses, generation of precursor metabolites and energy and those involved in carbohydrate metabolic process. Conclusions These EST-SSR markers presented in this work can be efficiently used for genetic mapping studies of segregating sugarcane populations. The high Polymorphism Information Content (PIC) and Discriminant Power (DP) presented facilitate the QTL identification and marker-assisted selection due the association with functional regions of the genome became an important tool for the sugarcane breeding program.

Comparison of microsatellites and isozymes in genetic diversity studies of Oryza glumaepatula (Poaceae) populations

The study of the genetic structure of wild plant populations is essential for their management and conservation. Several DNA markers have been used in such studies, as well as isozyme markers. In order to provide a better comprehension of the results obtained and a comparison between markers which will help choose tools for future studies in natural populations of Oryza glumaepatula, a predominantly autogamous species, this study used both isozymes and microsatellites to assess the genetic diversity and genetic structure of 13 populations, pointing to similarities and divergences of each marker, and evaluating the relative importance of the results for studies of population genetics and conservation. A bulk sample for each population was obtained, by sampling two to three seeds of each plant, up to a set of 50 seeds. Amplified products of eight SSR loci were electrophoresed on non-denaturing polyacrylamide gels, and the fragments were visualized using silver staining procedure. Isozyme analyses were conducted in polyacrylamide gels, under a discontinuous system, using six enzymatic loci. SSR loci showed higher mean levels of genetic diversity (A=2.83, p=0.71, A(P)=3.17, H-o=0.081, H-e=0.351) than isozyme loci (A=1.20, p=0.20, A(P)=1.38, H-o=0.006, H-e=0.056). Interpopulation genetic differentiation detected by SSR loci (R-ST=0.631, equivalent to F-ST=0.533) was lower than that obtained with isozymes (F-ST=0.772). However, both markers showed high deviation from Hardy-Weinberg expectations (F-IS=0.744 and 0.899, respectively for SSR and isozymes). The mean apparent outcrossing rate for SSR ((t) over bar (a)=0.14) was higher than that obtained using isozymes ((t) over bar (a)=0.043), although both markers detected lower levels of outcrossing in Amazonia compared to the Pantanal. The migrant number estimation was also higher for SSR (Nm=0.219) than isozymes (Nm=0.074), although a small number for both markers was expected due to the mode of reproduction of this species, defined as mixed with predominance of self fertilization. No correlation was obtained between genetic and geographic distances with SSR, but a positive correlation was found between genetic and geographic distances with isozymes. We conclude that these markers are divergent in detecting genetic diversity parameters in O. glumaepatula and that microsatellites are powerful for detecting information at the intra-population level, while isozymes are more powerful for inter-population diversity, since clustering of populations agreed with the expectations based on the geographic distribution of the populations using this marker. Rev. Biol. Trop. 60 (4): 1463-1478. Epub 2012 December 01.

Genomic analysis of ERVWE2 locus in patients with multiple sclerosis: absence of genetic association but potential role of human endogenous retrovirus type W elements in molecular mimicry with myelin antigen

Human endogenous retroviruses (HERVs) arise from ancient infections of the host germline cells by exogenous retroviruses, constituting 8% of the human genome. Elevated level of envelope transcripts from HERVs-W has been detected in CSF, plasma and brain tissues from patients with Multiple Sclerosis (MS), most of them from Xq22.3, 15q21.3, and 6q21 chromosomes. However, since the locus Xq22.3 (ERVWE2) lack the 5' LTR promoter and the putative protein should be truncated due to a stop codon, we investigated the ERVWE2 genomic loci from 84 individuals, including MS patients with active HERV-W expression detected in PBMC. In addition, an automated search for promoter sequences in 20 kb nearby region of ERVWE2 reference sequence was performed. Several putative binding sites for cellular cofactors and enhancers were found, suggesting that transcription may occur via alternative promoters. However, ERVWE2 DNA sequencing of MS and healthy individuals revealed that all of them harbor a stop codon at site 39, undermining the expression of a full-length protein. Finally, since plaque formation in central nervous system (CNS) of MS patients is attributed to immunological mechanisms triggered by autoimmune attack against myelin, we also investigated the level of similarity between envelope protein and myelin oligodendrocyte glycoprotein (MOG). Comparison of the MOG to the envelope identified five retroviral regions similar to the Ig-like domain of MOG. Interestingly, one of them includes T and B cell epitopes, capable to induce T effector functions and circulating Abs in rats. In sum, although no DNA substitutions that would link ERVWE2 to the MS pathogeny was found, the similarity between the envelope protein to MOG extends the idea that ERVEW2 may be involved on the immunopathogenesis of MS, maybe facilitating the MOG recognizing by the immune system. Although awaiting experimental evidences, the data presented here may expand the scope of the endogenous retroviruses involvement on MS pathogenesis

Brazilian urban population genetic structure reveals a high degree of admixture

Advances in genotyping technologies have contributed to a better understanding of human population genetic structure and improved the analysis of association studies. To analyze patterns of human genetic variation in Brazil, we used SNP data from 1129 individuals - 138 from the urban population of Sao Paulo, Brazil, and 991 from 11 populations of the HapMap Project. Principal components analysis was performed on the SNPs common to these populations, to identify the composition and the number of SNPs needed to capture the genetic variation of them. Both admixture and local ancestry inference were performed in individuals of the Brazilian sample. Individuals from the Brazilian sample fell between Europeans, Mexicans, and Africans. Brazilians are suggested to have the highest internal genetic variation of sampled populations. Our results indicate, as expected, that the Brazilian sample analyzed descend from Amerindians, African, and/or European ancestors, but intermarriage between individuals of different ethnic origin had an important role in generating the broad genetic variation observed in the present-day population. The data support the notion that the Brazilian population, due to its high degree of admixture, can provide a valuable resource for strategies aiming at using admixture as a tool for mapping complex traits in humans. European Journal of Human Genetics (2012) 20, 111-116; doi:10.1038/ejhg.2011.144; published online 24 August 2011

ALDH1A2 (RALDH2) genetic variation in human congenital heart disease

Abstract Background Signaling by the vitamin A-derived morphogen retinoic acid (RA) is required at multiple steps of cardiac development. Since conversion of retinaldehyde to RA by retinaldehyde dehydrogenase type II (ALDH1A2, a.k.a RALDH2) is critical for cardiac development, we screened patients with congenital heart disease (CHDs) for genetic variation at the ALDH1A2 locus. Methods One-hundred and thirty-three CHD patients were screened for genetic variation at the ALDH1A2 locus through bi-directional sequencing. In addition, six SNPs (rs2704188, rs1441815, rs3784259, rs1530293, rs1899430) at the same locus were studied using a TDT-based association approach in 101 CHD trios. Observed mutations were modeled through molecular mechanics (MM) simulations using the AMBER 9 package, Sander and Pmemd programs. Sequence conservation of observed mutations was evaluated through phylogenetic tree construction from ungapped alignments containing ALDH8 s, ALDH1Ls, ALDH1 s and ALDH2 s. Trees were generated by the Neighbor Joining method. Variations potentially affecting splicing mechanisms were cloned and functional assays were designed to test splicing alterations using the pSPL3 splicing assay. Results We describe in Tetralogy of Fallot (TOF) the mutations Ala151Ser and Ile157Thr that change non-polar to polar residues at exon 4. Exon 4 encodes part of the highly-conserved tetramerization domain, a structural motif required for ALDH oligomerization. Molecular mechanics simulation studies of the two mutations indicate that they hinder tetramerization. We determined that the SNP rs16939660, previously associated with spina bifida and observed in patients with TOF, does not affect splicing. Moreover, association studies performed with classical models and with the transmission disequilibrium test (TDT) design using single marker genotype, or haplotype information do not show differences between cases and controls. Conclusion In summary, our screen indicates that ALDH1A2 genetic variation is present in TOF patients, suggesting a possible causal role for this gene in rare cases of human CHD, but does not support the hypothesis that variation at the ALDH1A2 locus is a significant modifier of the risk for CHD in humans.

The humankind genome: from genetic diversity to the origin of human diseases

Genome-wide association studies have failed to establish common variant risk for the majority of common human diseases. The underlying reasons for this failure are explained by recent studies of resequencing and comparison of over 1200 human genomes and 10 000 exomes, together with the delineation of DNA methylation patterns (epigenome) and full characterization of coding and noncoding RNAs (transcriptome) being transcribed. These studies have provided the most comprehensive catalogues of functional elements and genetic variants that are now available for global integrative analysis and experimental validation in prospective cohort studies. With these datasets, researchers will have unparalleled opportunities for the alignment, mining, and testing of hypotheses for the roles of specific genetic variants, including copy number variations, single nucleotide polymorphisms, and indels as the cause of specific phenotypes and diseases. Through the use of next-generation sequencing technologies for genotyping and standardized ontological annotation to systematically analyze the effects of genomic variation on humans and model organism phenotypes, we will be able to find candidate genes and new clues for disease’s etiology and treatment. This article describes essential concepts in genetics and genomic technologies as well as the emerging computational framework to comprehensively search websites and platforms available for the analysis and interpretation of genomic data.

Testing allele homogeneity: the problem of nested hypotheses

Background: The evaluation of associations between genotypes and diseases in a case-control framework plays an important role in genetic epidemiology. This paper focuses on the evaluation of the homogeneity of both genotypic and allelic frequencies. The traditional test that is used to check allelic homogeneity is known to be valid only under Hardy-Weinberg equilibrium, a property that may not hold in practice. Results: We first describe the flaws of the traditional (chi-squared) tests for both allelic and genotypic homogeneity. Besides the known problem of the allelic procedure, we show that whenever these tests are used, an incoherence may arise: sometimes the genotypic homogeneity hypothesis is not rejected, but the allelic hypothesis is. As we argue, this is logically impossible. Some methods that were recently proposed implicitly rely on the idea that this does not happen. In an attempt to correct this incoherence, we describe an alternative frequentist approach that is appropriate even when Hardy-Weinberg equilibrium does not hold. It is then shown that the problem remains and is intrinsic of frequentist procedures. Finally, we introduce the Full Bayesian Significance Test to test both hypotheses and prove that the incoherence cannot happen with these new tests. To illustrate this, all five tests are applied to real and simulated datasets. Using the celebrated power analysis, we show that the Bayesian method is comparable to the frequentist one and has the advantage of being coherent. Conclusions: Contrary to more traditional approaches, the Full Bayesian Significance Test for association studies provides a simple, coherent and powerful tool for detecting associations.

Differentiation of African Components of Ancestry to Stratify Groups in a Case-Control Study of a Brazilian Urban Population

Background: Balancing the subject composition of case and control groups to create homogenous ancestries between each group is essential for medical association studies. Methods: We explored the applicability of single-tube 34-plex ancestry informative markers (AIM) single nucleotide polymorphisms (SNPs) to estimate the African Component of Ancestry (ACA) to design a future case-control association study of a Brazilian urban sample. Results: One hundred eighty individuals (107 case group; 73 control group) self-described as white, brown-intermediate or black were selected. The proportions of the relative contribution of a variable number of ancestral population components were similar between case and control groups. Moreover, the case and control groups demonstrated similar distributions for ACA <0.25 and >0.50 categories. Notably a high number of outlier values (23 samples) were observed among individuals with ACA <0.25. These individuals presented a high probability of Native American and East Asian ancestral components; however, no individuals originally giving these self-described ancestries were observed in this study. Conclusions: The strategy proposed for the assessment of ancestry and adjustment of case and control groups for an association study is an important step for the proper construction of the study, particularly when subjects are taken from a complex urban population. This can be achieved using a straight forward multiplexed AIM-SNPs assay of highly discriminatory ancestry markers.

Association between the C242T polymorphism in the p22phox gene with arterial stiffness in the Brazilian population

de Oliveira Alvim R, Lima Santos PCJ, Goncalves Dias R, Rodrigues MV, de Sa Cunha R, Mill JG, Junior WN, Krieger JE, Pereira AC. Association between the C242T polymorphism in the p22phox gene with arterial stiffness in the Brazilian population. Physiol Genomics 44: 587-592, 2012. First published April 10, 2012; doi:10.1152/physiolgenomics.00122.2011.-NADPH oxidase p22phox subunit is responsible for the production of reactive oxygen species in the vascular tissue. The C242T polymorphism in the p22phox gene has been associated with diverse coronary artery disease phenotypes, but the findings about the protective or harmful effects of the T allele are still controversial. Our main aim was to assess the effect of p22phox C242T genotypes on arterial stiffness, a predictor of late morbidity and mortality, in individuals from the general population. We randomly selected 1,178 individuals from the general population of Vitoria City, Brazil. Genotypes for the C242T polymorphism were detected by PCR-RFLP, and pulse wave velocity (PWV) values were measured with a noninvasive automatic device Complior. p22phox and TNF-alpha gene expression were quantified by real-time PCR in human arterial mammary smooth muscle cells. In both the entire and nonhypertensive groups: individuals carrying the TT genotype had higher PWV values and higher risk for increased arterial stiffness [odds ratio (OR) 1.93, 95% confidence interval (CI) 1.27-2.92 and OR 1.78, 95% CI 1.07-2.95, respectively] compared with individuals carrying CC + CT genotypes, even after adjustment for covariates. No difference in the p22phox gene expression according C242T genotypes was observed. However, TNF-alpha gene expression was higher in cells from individual carrying the T allele, suggesting that this genetic marker is associated with functional phenotypes at the gene expression level. In conclusion, we suggest that p22phox C242T polymorphism is associated with arterial stiffness evaluated by PWV in the general population. This genetic association shed light on the understanding of the genetic modulation on vascular dysfunction mediated by NADPH oxidase.

MYLIP p.N342S polymorphism is not associated with lipid profile in the Brazilian population

Background: A recent study investigated the MYLIP region in the Mexican population in order to fine-map the actual susceptibility variants of this locus. The p.N342S polymorphism was identified as the underlying functional variant accounting for one of the previous signals of genome-wide association studies and the N342 allele was associated with higher cholesterol concentrations in Mexican dyslipidemic individuals. To date, there is no further evaluation on this genotype-phenotype association in the literature. In this scenario, and because of a possible pharmacotherapeutic target of dyslipidemia, the main aim of this study was to assess the influence of the MYLIP p.N342S polymorphism on lipid profile in Brazilian individuals. Methods: 1295 subjects of the general population and 1425 consecutive patients submitted to coronary angiography were selected. General characteristics, biochemical tests, blood pressures, pulse wave velocity, and coronary artery disease scores were analyzed. Genotypes for the MYLIP rs9370867 (p.N342S, c.G1025A) polymorphism were detected by high resolution melting analysis. Results: No association of the MYLIP rs9370867 genotypes with lipid profile, hemodynamic data, and coronary angiographic data was found. Analysis stratified by hyperlipidemia, gender, and ethnicity was also performed and the sub-groups presented similar results. In both general population and patient samples, the MYLIP rs9370867 polymorphism was differently distributed according to ethnicity. In the general population, subjects carrying GG genotypes had higher systolic blood pressure (BP), diastolic BP, and mean BP values (129.0 +/- 23.3; 84.9 +/- 14.6; 99.5 +/- 16.8 mmHg) compared with subjects carrying AA genotypes (123.7 +/- 19.5; 81.6 +/- 11.8; 95.6 +/- 13.6 mmHg) (p = 0.01; p = 0.02; p = 0.01, respectively), even after adjustment for covariates. However, in analysis stratified by ethnicity, this finding was not found and there is no evidence that the polymorphism influences BP. Conclusion: Our findings indicate that association studies involving this MYLIP variant can present distinct results according to the studied population. In this moment, further studies are needed to reaffirm if the MYLIP p.N342S polymorphism is functional or not, and to identify other functional markers within this gene.

Using Prior Information from the Medical Literature in GWAS of Oral Cancer Identifies Novel Susceptibility Variant on Chromosome 4-the AdAPT Method

Background: Genome-wide association studies (GWAS) require large sample sizes to obtain adequate statistical power, but it may be possible to increase the power by incorporating complementary data. In this study we investigated the feasibility of automatically retrieving information from the medical literature and leveraging this information in GWAS. Methods: We developed a method that searches through PubMed abstracts for pre-assigned keywords and key concepts, and uses this information to assign prior probabilities of association for each single nucleotide polymorphism (SNP) with the phenotype of interest - the Adjusting Association Priors with Text (AdAPT) method. Association results from a GWAS can subsequently be ranked in the context of these priors using the Bayes False Discovery Probability (BFDP) framework. We initially tested AdAPT by comparing rankings of known susceptibility alleles in a previous lung cancer GWAS, and subsequently applied it in a two-phase GWAS of oral cancer. Results: Known lung cancer susceptibility SNPs were consistently ranked higher by AdAPT BFDPs than by p-values. In the oral cancer GWAS, we sought to replicate the top five SNPs as ranked by AdAPT BFDPs, of which rs991316, located in the ADH gene region of 4q23, displayed a statistically significant association with oral cancer risk in the replication phase (per-rare-allele log additive p-value [p(trend)] = 2.5 x 10(-3)). The combined OR for having one additional rare allele was 0.83 (95% CI: 0.76-0.90), and this association was independent of previously identified susceptibility SNPs that are associated with overall UADT cancer in this gene region. We also investigated if rs991316 was associated with other cancers of the upper aerodigestive tract (UADT), but no additional association signal was found. Conclusion: This study highlights the potential utility of systematically incorporating prior knowledge from the medical literature in genome-wide analyses using the AdAPT methodology. AdAPT is available online (url: http://services.gate.ac.uk/lld/gwas/service/config).

Relationship between Vitamin D Receptor gene polymorphisms and the components of metabolic syndrome

Abstract Background The Vitamin D Receptor gene (VDR) is expressed in many tissues and modulates the expression of several other genes. The purpose of this study was to investigate the association between metabolic syndrome (MetSyn) with the presence of VDR 2228570 C > T and VDR 1544410 A > G polymorphisms in Brazilian adults. Methods Two hundred forty three (243) individuals were included in a cross-sectional study. MetSyn was classified using the criteria proposed by National Cholesterol Educational Program - Adult Treatment Panel III. Insulin resistance and β cell secretion were estimated by the mathematical models of HOMA IR and β, respectively. The VDR 2228570 C > T and VDR 1544410 A > G polymorphisms were detected by enzymatic digestion and confirmed by allele specific PCR or amplification of refractory mutation. Results Individuals with MetSyn and heterozygosis for VDR 2228570 C > T have higher concentrations of iPTH and HOMA β than those without this polymorphism, and subjects with recessive homozygosis for the same polymorphisms presented higher insulin resistance than those with the heterozygous genotype. There is no association among VDR 1544410 A > G and components of MetSyn, HOMA IR and β, serum vitamin D (25(OH)D3) and intact parathormone (iPTH) levels in patients with MetSyn. A significant lower concentration of 25(OH)D3 was observed only in individuals without MetSyn in the VDR 1544410 A > G genotype. Additionally, individuals without MetSyn and heterozygosis for VDR 2228570 C > T presented higher concentration of triglycerides and lower HDL than those without this polymorphism. Conclusions Using two common VDR polymorphism data suggests they may influence insulin secretion, insulin resistance an serum HDL-cholesterol in our highly heterogeneous population. Whether VDR polymorphism may influence the severity of MetSyn component disorder, warrants examination in larger cohorts used for genome-wide association studies.